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^^ Abstract. The COMPASS experiment at CERN is designed for light hadron spectroscopy with emphasis on the detection of 

^"^ new states, in particular the search for exotic states and glue-balls. After a short pilot run in 2004 (190 GeV/c negative pion 

^^ beam, lead target) showing significant production strength for an exotic J^^ = 1~+ state at 1.66GeV/c^, we have collected 

Cn data with a 190 GeV/c negative charged hadron beam on a proton (liquid hydrogen) and nuclear targets in 2008 and 2009. 

1— H The spectrometer features good coverage by electromagnetic calorimetry, and our data provide excellent opportunity for 

^ simultaneous observation of new states in two different decay modes in the same experiment. The diffractively produced 

^^ (37r)~ system for example can be studied in both modes 7i~ p -^ 7i~7i^7i~ p and 7i~ p ^ 7i~7fi7fi p. Charged and neutral 

[-^ mode rely on completely different parts of the spectrometer. Observing a new state in both modes provides important cross- 

^-H check. First results of a preliminary PWA performed on the 2008 data are presented. 

^^^ Keywords: Hadron spectroscopy, light meson spectrum, gluonic excitations, exotic mesons, hybrids 

'T7' PACS: 13.25.-k,13.85.-t,14.40.Be,29.30.-h 

^ INTRODUCTION 

I— I The COMPASS fixed target experiment [1] at CERN SPS is dedicated to the study of nucleon spin structure and hadron 

spectroscopy, addressing the question of how nucleons and hadrons in general are built up from quarks and gluons [2] . 

K^ The COMPASS Collaboration has already collected data scattering a polarised 160GeV/c muon beam on polarised 

^ deuteron (^LiD) and proton (NH3) targets during the years 2002-2004 and 2006-2007. The gluon contribution to the 

\f^ nucleon spin is one example of physics determined from these data. During a second phase dedicated to physics with 

^\ hadron beams, we have collected unprecedented statistics of data with 190 GeV/c charged hadron beams on a proton 

CN and nuclear targets in 2008 and 2009. The feasibility of our apparatus for light mesons spectroscopy has been studied 

f*"^ in a short pilot run in 2004 (190 GeV/c negative pion beam, lead target). Based on the few days diffractive pion data 

^^ included in the 2004 run, pion dissociation into TC~ 7l~ 71^ final states has been analysed showing significant production 

^^ strength for an exotic J^^ = 1~+ state at 1.66 GeV/c^, which can be interpreted as the ;ri (1600) [3]. 

^"^ The high statistics data sample taken with the improved spectrometer in 2008/09 allows us not only to complete 

^ the search for the Tli but also to extent our analyses to further channels of interest (with lower cross section and 

V^ higher masses) and further develop our PWA methods. In particular the detection of final states with both charged and 

rN neutral particles is one of the key advantages of COMPASS as compared to previous fixed target experiments. First 

^ preliminary results on the 2008 data for pion dissociation into 3 pion final states, neutral mode: ;r~ n^n^, are presented. 

The simultaneous study of both modes allows for important cross-check (acceptances, systematics) and independent 
confirmation of any new state observed in the charged mode. 

LIGHT MESON SPECTROSCOPY - DIFFRACTIVE DISSOCIATION 

The naive Constituent Quark Model (CQM) characterises mesons as bound colour- singlet states of a quark q and a 
anti-quark q with flavours u^d and s grouped into SU(3)fiavour multiplets. Their total angular momentum /, parity P 
and charge conjugation C are given by / = \L — S\...\L-\-S\, P = (— 1)^+^ and C = (—1)^+*^ , where L is the relative 
orbital angular momentum of q and q, and S the total intrinsic spin (5" = 0, 1) of the qq pair. In addition the isospin / 
and the G parity defined as G = (—1)^+^+*^ are introduced, also conserved in strong interactions. Given the simplicity, 
the CQM is astonishingly successful in describing part of the meson properties as well as - to a large extend - the 
observed spectrum. 
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FIGURE 1. Left: (top) Meson production in diffractive scattering via t-channel Reggeon exchange, (bottom) Diffractive dissoci- 
ation into 371 final states as described in the isobar model: The diffractively produced resonance X~ with quantum numbers J^^M^ 
decays into an isobar with spin S and relative orbital angular momentum L with respect to the TTbacheior^ the isobar subsequently 
decays into two pions. At high energies, the Pomeron is the dominant Regge-trajectory. Right: Sketch of the two-stage COMPASS 
spectrometer (^ 50 m long) as used during hadron runs 2008 and 2009. 



In QCD, however, interactions between coloured quarks are described by exchange of gluons g carrying colour 
themselves, resulting in the prediction of new phenomena. In particular, colour- singlet mesons are not restricted to 
be composed of qq pairs but may consist of other colour-neutral configurations, like e.g. qqqq (tetra-quarks), qqg 
(hybrids) or gg (glueballs), which are mostly discussed in this context. Due to mixing with ordinary qq states with 
same quantum numbers /^^, such configurations are extremely difficult to find experimentally, since it is hardly 
possible to disentangle the contribution of each configuration. The experimental observation of spin-exotic mesons 
with quantum numbers forbidden in the CQM, like e.g. J^^ = , 0^~ , 1 ^, would thus provide a clear evidence for 
physics beyond the naive quark model and a fundamental confirmation of QCD. 

The lowest-lying hybrid is expected to have J^^ = 1~+. Lattice-QCD simulations [4] and flux-tube model calcu- 
lations [5] predict a mass between 1.7 and 2.2GeV/c^, and a preferred decay into bin and fiTi. Experimentally two 
candidates for a 1~+ hybrid have been found, 7ri(1400) and 7ri(1600), however, both are still heavily disputed in 
the community. The ;ri(1400) was mainly seen in rjTi decays, by e.g. E852 [6], VES[7], and Crystal Barrel [8]. The 
7ri(1600) was observed by both E852 and VES in the decay channels: pn [9, 10], ri^n [7, 11], fiTi [12, 13], and 
COJITT [13, 14]. Especially the observations of 7ri(1600) into pn based on analyses of 7i~7r^7i~ final state events are 
controversially discussed [15, 16]. 

In diffractive pion dissociation (at high energy), see Fig. 1 (left/top), the incident beam particle a is excited via 
(t-channel) Reggeon exchange to some resonance c, which further dissociates into n final state particles, whereas the 
target particle b remains intact: a + Z? ^ c + J, with c ^ 1 + ... + n particles, and d denotes the recoil (target) particle. 
Interactions of this type are characterised by two kinematic variables s and ^' = |^| — |^|min. where s = {p^^ p\^Y \^ 
the squared centre-of-mass energy and t = {p^, — PcY i^ the square of the four momentum transfered from the incident 
beam particle to the outgoing system c. Depending on the produced invariant mass mc, a minimum value of |^|niin is 
allowed by kinematics, which is small but larger than zero due to the longitudinal four-momentum transfer needed 
(mc > ma). In the centre-of-mass system: 



t' = t 



k|min = 2|^a||^c|(l-COSeo) >0 with \t {min = 2{E^Ec - \pa\\Pb\) - {^l^ ml) 



(1) 



where Oq is the scattering angle. Diffractive reactions have a total cross section in the order of 1-2 mb. Even though 
the differential cross section drops as 1/mc, states beyond 3GeV/c^ can be produced diffractively in a fixed target 
experiment like COMPASS (190 GeV/c^ 7l beam, proton target). Due to the forward kinematics the final state particles 
have to be detected mostly under small angles (with respect to the beam) requiring excellent angular resolution. 
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FIGURE 2. Exclusive events are selected by three mains cuts: A<^, exclusivity, and n^ mass. Left: A<^ vs. calculated beam energy 
(i.e. exclusivity). Centre: Exclusivity before/after A<^ cut. Right: Invariant mass of 7172 vs. 7374, cf. discussion in text. 

EXPERIMENTAL SET-UP IN 2008/09 

A detailed description of the COMPASS two-stage spectrometer (Fig. 1 (right)) dedicated to a variety of fixed- target 
physics programmes can be found in [17]. For the measurement with hadron beams started in 2008, a 40 cm long 
liquid hydrogen target with a diameter of 35 mm, or simple disks of solid material (part of 2009 run) have been used. 
The spectrometer features electromagnetic and hadronic calorimetry in both stages. Photon detection in a wide angular 
range with high resolution is crucial for decay channels involving tt^, T] or T]'. Therefore, the read-out electronics have 
been upgraded (from 10 to 12 bit SADCs) in 2008, allowing for Digital Signal Processing of the ADC signals, and 800 
of the lead glass Cherenkov counters (3000 in total) have been replaced by so-called Shashlik sampling calorimeters 
in the central part of ECAL2 to improve the radiation hardness as well as the energy resolution, see e.g. [18]. A new 
monitoring laser system for improved gain control of ECALl has further been installed in 2009. A Recoil Proton 
Detector (RPD) consisting of 2 concentric barrels of scintillator slats read out by PMTs was introduced to trigger on 
interactions inside the target and to detect the recoil particle. It performs a time-of-flight measurement at high accuracy 
(^ 350 ps). Finally, two CEDAR were installed, to separate the kaons (~ 2.5 %) in the beam from the pions (or, in 
case of proton beam, the pion contribution from protons). 

DATA SELECTION 

The data presently analysed for diffractively produced 7t~7t^TC^ final states corresponds to ~ 10 % of the data taken 
with pion beam in 2008. The diffractive trigger selected events with one incoming charged beam particle and a recoil 
proton detected by the RPD. Non-interacting beam and events out of acceptance were vetoed. Exactly one primary 
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FIGURE 3. Invariant mass spectra of - Left: Total outgoing system. Centre: 7Z it system. Right: tz^tt system 
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FIGURE 4. Left: Distribution of squared four-momentum transfer t' (Eq.l). Centre: Dalitz plot in ^2 region (1.320 ± 
O.lOOGeV/c^). Right: Dalitz plot in 712 region (1.670 ± 0.100 GeV/c^). 

vertex inside the target volume is required for each event. Events with exactly one outgoing charged track and 4 7s, 
from the two 71^ decays, detected in ECALl and ECAL2 were selected, if they give exactly one 71^ 71^ combination 
within a circular cut of ± 20MeV/c^ around the PDG mass (preselection, later tightened). Background events from 
elastic scattering have been suppressed by cutting on the energy Ej^.- < 185 GeV. In order to select exclusive events, 
three main cuts are applied consistently in terms of ±2 a of each distribution of the three observables after having 
applied the other two. Those are the angle A<l> (±0.2 rad), defined as the azimuthal angle between the total momentum 
of the outgoing pion system measured with the spectrometer and the one of the recoil proton detected with the RPD, 
which should be anti-parallel by momentum conservation, the exclusivity (±6 GeV) applied on the beam energy 
calculated from the outgoing system under assumption of energy conservation, and the 71^ mass (±16MeV). The 
corresponding distributions are shown in Fig. 2 (left), the exclusive sample appears around the nominal beam energy. 
The background below the exclusivity peak is suppressed by applying the cut on A<l> as clearly seen in Fig. 2 (centre). 
The final jiTi versus 7374 invariant mass distribution after all three 2 a cuts is given by Fig. 2 (right). It should be 
noticed that demanding exactly 4 ECAL clusters reduces presently the statistics outcome significantly, however, this 
will improve once our electromagnetic reconstruction is finalised, taking full advantage of the detector upgrade. The 
resultant invariant mass distributions are given in Fig. 3: The total 371 system (left) looks similar to the charged mode 
[19], clearly visible are the prominent ai(1260), ^2(1320) and 7r2(1670). The p~ is cleanly seen in the 7r~7r^ mass 
spectrum (centre) as well as the /2(1270) in 7r~7r^ (right), also we might see the /o(980) in the 71^71^ spectrum. Fig. 4 
(left) shows the t' distribution for the final sample of 240 k events. In Fig.4 (centre/right) the Dalitz plots in the ^2 
and 712 region, respectively, are shown. The main decays into the p~(770) are seen including the effect of constructive 
interference. Otherwise for the 712 region, the /2(1270) is not yet visible due to lack of statistics. 



PARTIAL WAVE ANALYSIS (MASS INDEPENDENT) 

A PWA has been performed restricted to the range 0.1 < ^' < 1.0, to stay above the RPD threshold and to ensure 
diffractive reactions. In order to determine all resonances present in the data, including the quantum numbers, we 
perform our PWA in two steps: a mass independent PWA and a subsequent mass dependent fit. The program used 
was originally developed in Illinois and modified at Protvino and Munich [20, 21]. At this first glance, essentially 
the same model that was used to analyse the 2004 data [3] is applied to the 2008 data to the neutral and charged 
[19] decay modes started to be analysed. Five isobars, the p(770), /o(980), /2(1270), p3(1690), and the {7i7i)s, as a 
parameterisation of the broad (7(600) and /o(1370) [22], are included. In total 42 partial waves are fitted to the data, 
including a background wave, which is flat in the relevant Gottfried- Jackson (GJ) angles and added incoherently. We 
discuss the first preliminary mass-independent fits, where the angular distributions are fitted in 40 MeV/c^ bins of the 
371 invariant mass m = mc, assuming the production strength for a given wave to be constant within a given mass bin. 
The underlying formalism is based on two assumptions. The total cross-section is separated into a resonance and a 
recoil vertex, and the isobar model, see Fig. 1 (left/bottom), is used to describe the X~ decay into three pions as a 
two-step process, without any further final state interactions among the pions nor with the target. The decay goes via 
an intermediate 7 = di-pion resonance, the so-called isobar, decaying into a pion pair, and a so-called bachelor pion. 
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FIGURE 5. Comparison of PWA intensities of main waves for neutral vs. charged mode. Left: Intensities of the ^2 (2^+1+ going 



The isobar spin S and relative orbital angular momentum L between the isobar and the bachelor pion couple to the 
spin / of the resonance X~ . The 371 system has isospin / > in general, and we can assume / = 1 as no flavour-exotic 
mesons are known in the light quark sector. Since the final state under study comprises an odd number of pions (and 
thus negative G-parity), the charge conjugation is positive. The amplitudes are constructed in the reflectivity basis 
[23] so that the X~ spin projection M > and the reflectivity £ = ±1 describes the symmetry under reflection at the 
production plane. Amplitudes of different reflectivities do not interfere due to parity conservation. Moreover, at high 
^/s the reflectivity equals naturality of the exchanged Reggeon, e = -\-l corresponds to natural parity exchange like 
e.g. Pomeron mediated reactions. The full set of quantum numbers J^^M^[isobar 7l]L defines a partial wave, whereas 
/ and G are not explicitly specified since they are fixed by the incoming pion to /^ = 1 ~ . For the spin density matrix, 
we allow for rank Nr=2 to account for helicity flip and non-flip amplitudes at the baryon vertex, assuming the target 
nucleon stays intact. Finally, the observed intensities are parameterised as a coherent and incoherent sum over the 
partial wave amplitudes [23]: 



(^indQpi^.Jn) 



Nr 

II 

e=±lr=l 



X7;.f^f(T,m)/j||^f(T',m)|W 



(2) 



where the three body kinematics are described completely by five phase space coordinates represented by T, measured 
for each event. They are the input for calculating the decay amplitudes, ^f for each partial wave /, using the D-function 
formalism in the helicity frame. The complex numbers T-f, the so-called production amplitudes, contain information 
of strength and interference of the waves. They are obtained using an extended maximum-likelihood method. The 
spectrometer acceptance can be taken into account directly in this procedure. 



FIRST RESULTS, CHECK OF ISOSPIN SYMMETRY 

As mentioned before, simultaneous observation of both (3;r)~ modes in the same experiment provides important 
cross-check. When an isospin 1 resonance X~ is produced, and subsequently decays via an intermediate 7 = di-pion 
resonance, the yield in the neutral mode should be half of that in the charged mode. Otherwise, if the di-pion is an 
isovector, equal yields are expected for both modes. This isospin symmetry holds only if the branchings are completely 
determined by the Clebsch-Gordan coefficients. In general, however, Bose-Symmetrisation with the bachelor pion is 
obligatory and might affect the observed branchings. We checked by calculation (using the wave functions) that the 
expectation for observed intensities as formulated before are correct for all isobars going to pn. Here, the effect 
might indeed not be negligible, but is the same for charged and neutral mode, and therefore cancels out. On the other 
hand, the isospin symmetry needs to be modified for isobars going into /o,2^ due to interference effects from Bose- 
Symmetrisation. For example, in case of 0~+ into /o( 1400) TT S wave, this effect doubles the expected suppression 
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Comparison of PWA intensities of main waves for neutral vs. charged mode. Left: (712) 2+0+ into /2(1270) n S 
(712) 2~+0+ into p~7r F wave. 



factor of two (simply expected from Clebsch-Gordan coefficients) of intensities observed in the neutral versus the 
charged mode, see Tab.l. In this extreme case, it is rather a factor of 4. Otherwise we find no distortion of the isospin 
symmetry due to such interference effects for the example of 2~+ into /2(1270);r S wave (at the ;r2(1670) PDG mass), 
and expect here indeed to observe the pure suppression factor of two given by the Clebsch-Gordan coefficients. The 
calculated branching ratios summarised in Tab. 1 are in good agreement with the data. Even though acceptance has not 
yet been corrected for at this stage of the analysis, the dominant intensities show the symmetry as described above, 
and thus the acceptance is proved to be rather uniform. In order to compare the observed intensities in the neutral to 
the charged mode, the PWA results have been normalised to the well-established narrow 02(1 320) observed in both 
modes, see Fig. 5 (left). Looking at ai, l++0+[p~(770);r] S wave (Fig. 5 (right)), we find the intensities as well as the 
widths being quite similar for the different modes, as expected. For the 712, 2~+0+[/2(1270);r] S wave (Fig. 6 (left)), 
we obtain the neutral mode being suppressed by a factor ^2.2 relatively to the charged case, which is qualitatively 
already in good agreement with our expectation. On the other hand looking at 712, 2~+0+[p~(770);r] F wave (Fig. 6 
(right)), we find consistently again about the same intensities. As a quality check of the fits, we compare real data 
to Monte Carlo (MC) events, which were generated by weighting phase space events with the production and decay 
amplitudes from the fit result under the assumption of a uniform acceptance. Such comparison for the decay angles 
(cos 6 and (j)) of the p in the GJ frame are depicted for the neutral mode data in Fig. 7 (left). The angles are shown for 
the limited mass region around ai and a2 (L22 to L38 GeV/c^). Comparing the angles to the corresponding ones of the 
charged mode data. Fig. 7 (right), limited to the same mass range around the (22, one finds similar angular distributions, 
which is expected, since the physics is the same. Furthermore, for both cases the assumption of a flat acceptance seems 
to be valid, and therefore the comparison of both modes in terms of isospin symmetry is reasonable, even though the 
data has not yet been corrected for acceptance. 

TABLE 1. Isospin symmetry checks: Calculation of branching ratios (BR) for 
charged and neutral mode and different isobar decays. 

BR = N(7r~7r^7r^ )/N(7r~7r~7r+ ) - calculated from isobar model amplitudes 



BR(p7r) = l. 

BR( 0-+/o(1400)7r S) = 0.26 (at 1.3 GeV) = 0.29 (at 1.8 GeV) 

BR( 0-+/o(980)7r S) = 0.44 (at 1.8 GeV) 

BR( l++/o(1400)7r P) = 0.80 (at 1.3 GeV) 

BR( 2-+/2(1270);r S) = 0.50 (at 1.67 GeV) 
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SUMMARY & OUTLOOK 

The COMPASS experiment has a high potential for contributing to light meson spectroscopy. Data with charged hadron 
beams on different targets with high statistics have been taken in 2008/09 with an upgraded apparatus. One main goal 
of the spectroscopy program is the search for J^^ exotic states with gluonic degree of freedom, and to illuminate e.g. 
the disputed hybrid candidate ;ri (1600). In particular the detection of final states comprising both charged and neutral 
particles, allowing for cross-check and independent confirmation of any new state found, makes COMPASS unique 
as compared to previous fixed-target experiments. A first event selection and partial wave analysis of diffractively 
produced Tr'n^n^ final states of a subset of the 2008 data (pion beam, proton target) has been performed. The 
observed main waves are, at this first glance, in good agreement with theoretical expectations (isospin symmetry), 
which demonstrates the feasibility of COMPASS for hadron spectroscopy not only of charged but also of channels 
involving neutral particles. The data recorded is of sufficient statistics to even study systematically the isobar model 
itself. Next steps for this analysis are increasing the statistics, application of acceptance corrections, extension of the 
waveset, and studying the existence of the exotic 1~+ wave in the 2008/09 data set. 
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